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Abstract. A thermo-mechanical model describing hydrogen storage by use of metal hydrides has 
been recently proposed in fj] describing the formation of hydrides using the phase transition approach. 
By virtue of the laws of continuum thcrmo-mechanics, the model leads to a phase transition problem in 
terms of three state variables: the temperature, the phase parameter representing the fraction of one 
solid phase, and the pressure, and is derived within a generalization of the principle of virtual powers 



proposed by Fremond accounting for micro- forces, responsible for the phase transition, in the 
\ whole energy balance of the system. Three coupled nonlinear partial differential equations combined 

' with initial and boundary conditions have to be solved. The main difficulty in investigating the 

resulting system of partial differential equations relies on the presence of the squared time derivative 
p ^ ' of the order parameter in the energy balance equation, and actually this term was neglected in the 

analysis performed in [3]. Here, the global existence of a solution to the full problem is proved 
by exploiting known and sharp estimates on parabolic equations with right hand side in L . Some 
complementary results on stability and steady state solutions are also given. 



1. Introduction 

(N 

^ ■ Let us consider a bounded smooth open set Q C R , 1 < N < 3, with boundary T, and some 
final time T > 0. In the space-time cylindrical domain Q := Q x (0,T) we deal with the following 
evolution problem 

2 : d t e-A6 = -h(6)d tX + ti(d t x) 2 , (1.1) 

e = ,j>{e,x):=e-x(h{e)-eh'(e)), (1.2) 

lidtx-&X + Z = h(6)- logu, (1.3) 

>< : C e /3(x) := 9I m ( X ) + log(l + X), (1-4) 

F? I d t u-Ap = 0, (1.5) 

p = u{l + X ), (1.6) 

with the boundary conditions 

d n e = d nX = d nP + 1P = ( 7 >0), (1.7) 

in E = T x (0,T), and the initial conditions 

e(-,0) = e , X (-,0) = Xo, u(-,0) = u in a (1.8) 

The system (ll.ip - fll.8p originates from a thermo-mechanical model introduced in [1] and describing 
hydrogen storage by use of metal hydrides. Hydrogen storage, which turns out a challenging subject 
in energetic and industrial applications, basically implies a reduction in the enormous volume of 
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hydrogen gas: so pressure plays a fundamental role in the evolution of the phenomenon. There are 
basically six methods to store hydrogen reversibly with a high volumetric and gravimetric density: 
one of them is by metal hydrides. This technique exploits the possibility of many metals to absorb 
hydrogen: such metals and alloys are able to react spontaneously with hydrogen and they can store 
a large amount of it. These materials, either a defined compound or a solid solution, are designed as 
metal hydrides: their use in the hydrogen storage is of interest in terms of safety, global yield, and 
long-time storage. 

Following a usual approach in thermodynamics, Bonetti, Fremond and Lexcellent [I] derived the 
governing equilibrium equations (actually simplifying some physical assumptions, but ensuring ther- 
modynamical consistency). The model has next been resumed and extended in [8]. 

Let us comment about the physical meaning and derivation of equations ( 11. lh - ( 1131) . As already 
mentioned, metals are able to absorb hydrogen atoms and combine with them to form solid solutions. 
The existence of two solid phases is supposed and the occurrence of one phase with respect to the 
other depends on the pressure p of the hydrogen. The volume fraction x of one of the phases is taken 
as a state quantity. Then, assuming that no voids appear in the mixture, the volume fraction of the 
other hydrogen phase is simply given by 1 — x- Hence, x nas to satisfy the relation < x < 1> which 
is ensured by the inclusion (ll.4p due to the presence of the subdifferential of the indicator function 
Z[ 0; i] of the interval [0, 1] as a component of the maximal monotone graph /3. Recall that the indicator 
function /[ 0jl ] is given by 7[ 0i i](x) = if x £ [0, 1] , and I[o,i](x) — +°° otherwise. Note that (3 also 
contains the contribution of the smooth (when restricted to [0, 1]) increasing function log(l + x) (i n 
our notation, logr = s if r = e s ). The graph /3 is actually multivalued and maximal monotone (see, 
e.g., [3 [1]), with £ G fl(x) if and only if 

{< if x = 
= if < x < 1 . 
> if x = 1 

Besides the phase parameter x along with its gradient Vx accounting for local interactions between 
the different phases, the other state variables of the model are the absolute temperature 9 and the 
hydrogen pressure p. Constitutive relations for the state quantities are chosen in such a way that 
the principles of thermodynamics are satisfied. 

The first equation ( II .ip is derived from the first principle of thermodynamics. It is obtained gener- 
alizing the principle of virtual powers as proposed by Fremond in [TT], i.e. including internal power of 
micro-forces (which are responsible for phase transitions) in the energy balance. In particular, note 
that in the right hand side, we have the quadratic nonlinearity {dtx) 2 standing for the internal power 
of dissipative micro-forces. The auxiliary variable e is prescribed in terms of 9 and x by relation 
( II. 2p via the smooth function ip. In fact, ip satisfies all the good properties we need (in particular, 
invertibility with respect to the first variable) since the function h is smooth and small enough (see 
[HE]): more precisely, we assume that there exists a constant Ch > 1 such that 

heC 2 (R)nW 2 '°°(R), 1 1 h 1 1 W 2,~ (M) + sup I rh'(r) | < c h (1.9) 

— < 1 + rh"(r)s < c h for all r G R and s e [0, 1]. (1.10) 

Equation (ll.3p . describing the evolution of the phase parameter, comes from the already mentioned 
generalization of the principle of virtual power in which micro-movements (related to the evolution of 
X as they are linked to changes of micro-structure) are included [11] . Equation (II. 5p is deduced from 
the usual continuity equation, in which the velocity is assumed to be proportional to the gradient 
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of the pressure. The structure of equation (11. 3p (along with (jl.4p ) is by now quite standard in the 
framework of phase field models: observe however that in the right hand side one has to handle 
the singular term logu (or alternatively logp). Then, one has to recover admissibility and enough 
regularity for logu from both fll.5D - fll.6jl and the initial datum no, that has to be strictly positive. 
The same can be repeated for eo if we want to guarantee that the absolute temperature 9 remains 
positive during the evolution, as it should be: but concerning (II. ip . we can be more flexible and 
develop our existence theory also in the general case, without sign restriction. 

About boundary conditions fll.7p . we explain that d n denotes the outward normal derivative on 
the boundary. Hence, Neumann homogeneous boundary conditions are set for 9 and Xi while the 
condition for p states that the hydrogen flux through the boundary is somehow proportional to 
the value of the interior pressure on T. Anyhow, in our contribution we treat the case 7 = too, 
corresponding to the case when the system is perfectly isolated. 

Now, let us briefly discuss the presence of dissipative quadratic nonlinearities and the coupling 
of phase transitions with pressure. Concerning phase changes with microscopic motions, there is a 
comprehensive literature on a class of models, some of them originating from the theory developed by 
Fremond [11]. In this respect, we first refer the reader to [6]: the resulting system of phase field type 
is characterized by the presence of new nonlinearities, like {d t x) 2 an d 9d t Xi which were not present 
in the classical formulation of phase field systems (not accounting for microscopic stresses). Several 
authors have dealt with this kind of problems and various situations have been analyzed. However, 
mainly for analytical difficulties due to the presence of nonlinearities, to our knowledge there is no 
global in time well-posedness result for the complete related system in the 3D (or 2D) case. A global 
existence result is proved in the ID setting [131 [15] or f° r a non-diffusive phase evolution [9]. Other 
results have been obtained for some regularized versions of the problem [3]. Weak solutions in terms 
of energy conservation (of Feireisl type) have been recently constructed in [10]. In addition, if we 
focus on the presence of the pressure, and thus on the continuity equation, note that some recent 
papers have dealt with phase transitions allowing different densities of the phases or adding the 
possibility of voids (see, for instance, [5J [T2j, where internal constraints on the pressure are given 
through indicator functions). 

In our system, the unpleasant nonlinearity in the right hand side of (II. ip includes p (d t x) 2 with 
p being a positive parameter. Note that the same p appears as the coefficient of d t x in (H-3p . Since 
our aim is showing the existence of a possibly weak solution to the initial boundary value problem 
fll.ip - fll.8p . the presence of p (d t x) 2 m (H-ip creates, as we have already pointed out, a difficulty from 
the mathematical point of view. For this reason, the authors of [I] suppressed it in the first analysis 
they carried out after the derivation of the model. That term is instead maintained in [8], where 
equation (II. 3p is replaced by 

pd tX - vAd t X ~ Ax + f = h(9) - logu, (1.11) 

where v stands for a small positive coefficient. The additional term —vAdtX actually contributes for 
a better regularity of dtX, and consequently the nonlinearity p (dtx) 2 to be handled in (II. ip turns 
out to be smoother. 

In the present contribution we can overcome the mentioned difficulty by exploiting an argument 
by Boccardo and Gallouet [2] for the study of second-order nonlinear parabolic equations with right 
hand side in L 1 . We prove the global- in-time existence of a weak solution to (jl.ip - fll.8p by taking 
(in Section [3]) an approximation of the problem based on (II. lip , with a small parameter v intended 
to go to zero in the limit procedure. Then we derive a priori estimates in Section H] and finally 
perform the passage to the limit: this is the subject of Section [5] which concludes the existence proof. 
Section [6] is devoted to the proof of some complementary results stated in Section |2] (along with the 
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main existence theorem) and regarding positivity of the variable 9, stability estimates in the case 
7 > and stationary solutions for the case 7 = 0. 



2. Main result 

In this section, we make precise assumptions on the data, specify our notion of solution to problem 
(11.11) — (11.80 and state our existence result. For convenience we set 

V := H\tt) H:=L 2 (tt), W := {v G H 2 ({1) : d n v = on T} , (2.1) 

and endow these spaces with their standard norms, for which we use the self-explanatory notations 
|| ■ \\h, II • || v, and || • \\w (the norm in L 2 (Q; R N ) being also denoted by || • \\h)- We remark that the 
embeddings W C V C H are compact, because Q is bounded and smooth. Since V is dense in H, 
we can identify H with a subspace of V in the usual way (i.e., so as to have that (u, v) = (u, v)h for 
every u G H and v G V); the embedding H C V is also compact. Let T > and set Q = Q x (0, T) 
and E = T x (0,T). 

We next recall that 

fi > 0, 7>0, h:R->R satisfies (jOjl -l OHll (2.2) 
and introduce the operators Ay, A : V — > V specified by 

(AjV,w) := / Vv-Vwdx + ^j vwdx Vv,wEV, A = A . (2.3) 
Jn Jr 

In view of the definition of ip : M x [0, 1] — > M given in ( II. 2p . it is straightforward to check that the 
partial derivatives d\ip, (with respect to first and second variable, respectively) satisfy 

— < diip(r, s) = 1 + rh"(r)s < c h , \d 2 ij)(r, s)\ < \h(r)\ + \rh!(r)\ < c h 

c h (2.4) 

for all r G R and s G [0, 1]. 

In particular, by looking at ( II. 2p again, we remark that for all e G R and x £ [0, 1] there exists a 
unique 9 G R fulfilling e = ip(9, x), and we denote by the inverse function of i/j with respect to 
the first variable. Clearly, the partial derivatives of satisfy 

- < d^-\r )S ) < c h , \d 2 ij-\r,s)\ <\d 1 ^\r,s)d 2 ^-\r,s),s)\ < {c h f 

c h (2.5) 

for all r£l and s G [0, 1]. 

In order to develop our existence theory, due to the presence of the term fi (dtx) 2 i n (U-ip we should 
give a meaning to (II. ip also when the right hand side is just in L X (Q). To this aim, let us extend 
the operator A in O from W l > q {VL) to {W l > q '(ft))' ', 



(Av, w) :- 



j Vv ■ Vwdx for all (v,w) G W hq {n) x W hq '(n), 

N + 2 1 1 ^ 

where 1 < q < , - H = 1, q' > N + 2. 

N + 1 q q' 



We now introduce the proper convex lower semicontinuous function : R — y [0, +00) defined by 



log(l + s)ds if < r < 1 
Jo ^ } ~ ~ , (2.7) 

+00 elsewhere 
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so that f3(r) > /3(0) = for r G M, and its sub differential is nothing but /3 (cf. f ll.4jl ). Hence, let us 
consider the functionals (cf. [HI formulas (58)-(66)]) 



f / (5(v)dx if we// and^MeL 1 ^) 
j H (t,) = ) J n , J( v ) = J H (u) if w e V, {2.1 

{ +00 if v e # but /3(f) £ L^fi) 



It is known that the inclusion w G /3(f) a.e. in f2 can be rewritten as w G OhJh{v) provided 
(v,w) G H x H or reinterpreted in terms of the sub differential dyyJ : V —> 2 V (different from 
the previous sub differential duJu '■ H — > 2 H ) as w G c^y/J^f) if (v,w) G V x V. Thus, this last 
inclusion can be seen as an extension of the previous ones whenever v G V. We note in particular 
that 

d H J H (v) = HH d vy ,J{y) for all v G V, (2.9) 

and announce that for our convenience we will use dyyJ in the variational formulation of f ll.ip - fll.8p 
to account for relation (JO 



Concerning equation (jl.5]l . we point out that we shall take advantage of a time- integrated version 
of it which is obtained with the help of the corresponding initial condition in f ll.8p . that is, 

u-A{l*p)=u , (2.10) 

where (l*p)(-,i)= / p( • , s)ds for t G (0, T). 

At this point, we can prescribe the conditions on the initial data. We assume that 

e G L 1 ^), xoey with /3( Xo ) G L 1 ^), ^ 
u G L 1 ^) n 1/' such that logu G L X (Q). 

Let us just notice that the requirement on xo exactly says that Xo £ D(J), i.e. Xo nes i n the domain 
of J, and this obviously entails < Xo < 1 a - e - m ^- The assumption on u implies that both u 
and log no ar e in L 1 (f2), whence in particular uq is strictly positive a.e. in Q. 

Here is our notion of weak solution to the problem fll.ip - fll.8p . 
Definition 2.1. We say that the sextuple (e,9,x,^,u,p) is a weak solution to problem f ll.ip - fll.8p if 

e,6E L q {0, T; W lt9 (Q)) for some q fulfilling Kq< j^^, (2.12) 

d t ee L 1 (0,T;(W 1 ' q '(fl))') with q' = — ^— > N + 2, (2.13) 

q- 1 

X eH\0,T;H)nL 2 (0,T;V), (2.14) 

£GL 2 (0,T;y'), (2.15) 

u G L°°(0, T; V) fl L 2 (0, T; H) and u > a.e. in Q, (2.16) 

pEL 2 (0,T;H) with 1 * p G L°°(0, T; V) (2.17) 
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and the following equations and conditions are satisfied: 

d t e + A9 = -h(9)dtx + Mx 2 in (W 1 '*^))', a.e. m (0,T), (2.18) 

e = i/>(6,x) a.e. inQ, (2.19) 

Hd t x + Ax + t, = h{6) -logu in V, a.e. in (0, T), (2.20) 

£(t) G fyW(x(t)) /or a.e. t G (0,T), (2.21) 

m + A 7 (1*p)=m in V, a.e. in (0,T), (2.22) 

p = «(1 + %) a.e. inQ, (2.23) 

e(0) = e in (W^'ity)', x(0) = Xo in H. (2.24) 

As a remark, note that the right hand side of (I2.18P is only in L X (Q), and we have L x (f2) C 
(W^'in))'. Indeed, seeing that q' > N, then W 1,q ' (Vt) is densely and compactly embedded in 
L°°(Q). We also point out that the inclusion (I2.2ip yields 

< x < 1 a.e. in Q, (2.25) 

as a consequence, since x(t) nas to belong to D(dvyJ) Q D(J) for a.e. t G (0, T). The positivity 
requirement on u G L 2 (Q) (cf. (I2.16P ) ensures that logw is well defined and measurable as a function: 
in addition, by comparing terms in (I2.20p we read that logw G L 2 (0,T;V) at least. The initial 
conditions in (I2.24p make sense due to (I2.12I) - (I2.13I) . which entail the continuity of e from [0,T] to 
(W 1 ' 11 '^))', and (EH}. 

It is time for us to state our existence result. 

Theorem 2.2. Under the assumptions (12. 2p and (12. lip . there exists a weak solution (e,6,x,£,,u,p) 
to problem (1 1 . 1 j) — ( TTTH|) in the sense of Definition \2.1[ Moreover, the further regularity properties 

e,6e L°°(0, T; L\n)) n L r (Q) n L«(0, T; W 1 *^)) 

/or a// r, o fulfilling 1 < r < 1 < g < j^^ , (2.26) 

XGC°([0,T];V0nL 2 (0,T;WO, (2.27) 

^i 2 (0,T;4 (2.28) 

u G L°°(0, T; L 1 (fi)) and log -a G L°°(0, T; /^(fi)) D L 2 (0, T; V) (2.29) 
aoW true for (e, x, £, p) • 

Remark 2.3. In view of f l2~9|) and (12721"]) . it turns out that (127281) ensures the validity of G 
9hJh{x{P)) f° r a - e - ^ e (0, , which in turn yields (11.41) . that is £ G f3(x) a.e. in Q. Thanks to 
(I2.27P and (12.291) , we infer from (12.201) that (11.31) is solved almost everywhere as well, which is not the 
case for the other equations (11.11) and (11.51) . In fact, for that we would need a additional regularity 
for (e, 9) and (u,p), respectively. Note however that equation (12.101) is satisfied a.e. in Q if u G H, 
since in this case 

1 *p G L 2 (0,T;H 2 (n)); (2.30) 

indeed, (I2.30p is a consequence of (I2.17p . A(l*p) G L 2 (0, T; H) and known elliptic regularity estimates 
which exploit the smoothness of the domain and of the boundary condition. Anyway, we would like 
to emphasize that (11.51) holds true in the sense of distributions in Q. 
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The next statement regards positivity of the variable 9, which is consistent with the fact that it 
represents the absolute temperature. 

Theorem 2.4. Assume that 

\og6oeL\tt), where 9 := V _1 (e , Xo), (2.31) 
in addition to (12.21) and (12.111) . Then the solution (e, 0, %, £, specified by Theorem \2.2\ satisfies 

\og6 G L°°(0, T; L 1 ^)) fl L 2 (0, T; V), ^Gi 2 (Q) (2.32) 

and in particular 9 > a.e. in Q. 

We also have two results concerning the large time behaviour in the case 7 > and the steady 
state solutions for the problem with 7 = 0. 

Theorem 2.5. Let 7 > and the assumptions of Theorem \2.2\ hold. Then, we have that 
u G L 2 (0, +00; L 2 (f2)) and ||w(t)||v decays to zero exponentially as t — > +00. 

In connection with this statement, note that u(t n + •) — > strongly in L 2 (0,T; H) for all T > 
and for any sequence t n — >■ +00. In some improper way, this kind of convergence to zero of u would 
entail that the right hand side of (12.201) tends somehow to +00, due to the presence of — logw. Thus, 
we expect that x(t) goes to 1, or becomes 1 from a certain time, which corresponds to one of the 
two solid hydrogen phases that becomes dominant at a long range. 



Theorem 2.6. Let 7 = 0. If (e,9,x,£,,u,p) is a stationary weak solution to (ll.ip - fll.8p . then the 
elements of the sextuple are all constant functions satisfying 

if h(9) — \ogp > then x = 1 an d p = 2u, 

if h(9) — logp = then < x ^ 1 an d P = (1 + X) u i 

if h(9) — logp < then x = and p = u. 

3. Approximation of the problem 

In order to prove Theorem 12.2} we introduce a regularizing term in (I2.20p . take advantage of the 
results shown by Chiodaroli in jH], prove uniform estimates and then pass to the limit with respect 
to the approximation parameter. 

We also regularize initial data and deal with a sequence of strictly positive values {7 n }neN to be 
used in the boundary condition for p (cf. (II. 7p ): more precisely, we assume that 

7„ > , 7„ — > 7, whence A ln — > A 1 in £(V; V'), as n — > 00. (3.1) 

Besides, we require the following conditions for the approximating initial data: 

Xon ^ W fl D(d v ,y J), and consequently (3(xon) £ L 1 ^), (3.2) 
e n, 9 0n e V, where 9 0n := ^ _1 (e n, Xon), (3.3) 

won, Pon G V, where p 0n := u 0n (l + xon), and logM 0n e L 1 (Q), — — G H. (3.4) 

POn 
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Clearly, e 0n , Xon, w 0n are smoother than the data e , Xo, u o satisfying assumption (12.111) . Morever, 
we also ask that 

e 0n -»■ e in L 1 ^), Xon Xo hi V, u 0n w in V as n-) oo, (3.5) 
llXonllw ^ Cn i H M on - logito»||Li(n) < c for a11 " ^ N. (3.6) 

Here and in the sequel, C denotes a positive constant, possibly depending on data but independent 
of n, that may change from line to line and even in the same chain of inequalities. Let us also note 
that we do not need to prescribe any boundedness property for {/3(xon)} in L 1 ^): indeed, this is 
ensured by ( 13. 21) . which implies < xon < 1 a.e. in Q and consequently (cf. (12.71) ) 



J(X0n) = \\P(X0n)\\ L i {u) < J \X0n\ log 2 dx < 

Moreover, the strong convergence of Xon to Xo stated in 1I3.5P also entails 

J(Xon) -> J(xo) as n -> oo. (3.7) 

Let us now briefly sketch how to construct the sequence of approximating initial conditions 
(eon, Xon, uon) satisfying f !3.2[) — A3. T[) . First, there is no problem in constructing a sequence {eo n } 
that fulfills ( 13 .3p and converges to Cq in L x (f2): it suffices to extend eo outside of Vt by zero and then 
regularize it by taking a convolution product over M, N with a sequence of mollifiers. 

As for {xon}, w e consider the sequence of solutions to 

Xon + -Axon = Xo in V . (3.8) 

n 

Since xo £ V, we have xon £ W and < xon < 1 a.e. in fl by the maximum principle. Moreover, we 
have that Xon € D(d V y>J): indeed, log(l + Xon) e 5y,^'^(Xon) as 



^(XOn) = / /3(X0n)cfe< / log(l + XOn) (XOn ~ Xo) dx + / P(Xo) dx 

Jn Jn Jn 

= (log(l + X0n),X0n - Xo) + J(Xo)- 

By the choice (13.81) . the first inequality in (13.61) holds true. To check it, it suffices to test (I3.8P by 
Xon first and then by Axon, integrating by parts in the right hand side. Furthermore, in singular 
perturbations like 1 13. 8ft a standard argument allows to deduce the strong convergence xon — > Xo by 
proving weak convergence plus convergence of norms. 

The same kind of approximation can be assumed for uq. In fact, let us take the solution uq h to 

u 0n + -Au 0n = u Q + - in V, (3.9) 

n n 

with the additional 1/n in the right hand side. It is straightforward to see that u 0n £ V and u 0n > 1/n 
a.e. in Q, the latter following from the comparison principle. Then 113 .4p follows and l/po n is even in 
V fl L°°(Q). Further, if we test 1 13.91) by (1 — 1/won) and observe that r4 (1 — 1/r) is the derivative 
of the positive and convex function r4r - logr, we find out that 



(U 0n ~ log Won) dx - f (u - log U ) dx < [ [1 — ) (u 0n - Uq) dx 

Jn Jn V u onJ 



< [ fl — J (U 0n - Uq) dx + - [ |VlogM 0n | 2 rfx < / - f 1 —)dx< — 

Jn\ u onJ n J n Jn n \ u onJ n 

as well as 

0< / (u 0n -logwon) dx < / (u -logu )dx+ 
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and the property ||M 0n — log it 0n Hz^m) < C turns out from (12. lip . 

The following proposition is a consequence of [8j Theorems 2.3 and 2.5] and it holds for all n > 1. 

Proposition 3.1. Under the assumptions (12. 2 I) and (13.21) - (13.41) . t/iere exists a sextuple 

(^n; ^n; Xn; £ru \] Pn) 

fulfilling 

e n , £ n G H\0,T;H)nL°°(0,T ] V), (3.10) 

X„ G ^ 1 ' oo (0,r;V)nL oo (0,T;^)nL oo (g), (3.11) 

e n 6i 2 (0,T;n, (3.12) 

w n G i? 1 (0, T; if) n L°°(0, T; and logw n G L°°(0, T; Z/^fi)) D L 2 (0, T; V), (3.13) 

p n G ZZ^O, T; if) n L°°(0, T; 1/) n L 2 (0, T; ZZ 2 (ft)) (3.14) 

and £/je following equations and conditions are satisfied: 

d t e n + Ad n = -h{d n )d tX n + ^{dtXn) 2 inV, a.e. in (0,T), (3.15) 
e n = ip(0n,Xn) a.e. in Q, (3.16) 

l^d t Xn + -A(d t Xn) + Axn + t,n = h{0 n ) - log u n in V', a.e. in (0, T), (3.17) 
n 

e„(t) G <9yyJ(XnW) M a.e. t G (0,T), (3.18) 

« n + A ln {l *p n ) = u 0n inV, a.e. in (0,T), (3.19) 

p n = u n {l + Xn) a.e. in Q, (3.20) 

e n (0) = e „ 2n 1/', x re (0) = xon. m H. (3.21) 

After the statement, we note that the regularity properties H3. lQj) — A3. 14[) entail (I2.12I) - (I2.17I) and 
(12T26|) - (12T271) . (l2~29l) . On the other hand, (l2~28l) is not granted by (l3~T2l : in fact, from (T37T7|) it just 
follows that 

n^AidtXn) + £n G L°°(0, T; if) + L 2 (0, T; V) 

(cf. (13.111) and (13.131) ). but A(d t Xn) is only in L°°(0,T; V'). An additional remark concerns equation 
(13TT9D which can be rewritten for u n and p n as 

d t u n + A 7n p n = in V, a.e. in (0, T) (3.22) 

along with the initial condition 

u n (0) = u 0n in H. (3.23) 

Note that ( 13. 22ft and ( I3.13l) -( l3.14l) imply that u n and p n solve ( 11 .5[) and the related boundary condi- 
tion almost everywhere. The same can be established for (13.151) which yields the analog of (11.11) a.e. 
in Q: indeed, (153U|) - (j3TTT]) entail A9 n G L 2 (0,T; H) so that 9 n G L 2 (0,T; W) too. As a last remark, 
we point out that (13.101) . (I3.16P and (11.91) allow us to rewrite (13 . 1 5[) as 



(1 + Xn0 n h"(e n ))d t e n + A9 n = -e n h'{9 n )d t Xn + KdtXnf mH, a.e. in (0,T). 



(3.24) 
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4. A PRIORI ESTIMATES 



Now, we carry out the uniform estimates, independent of n, and useful to pass to the limit. 

First a priori estimate. Test (13.171) by Xn and (I3.22p by (1 — l/u n ): this is formal but it can be 
made rigorous by using a truncation of the function r i— > (1 — 1/r) in a neighborhood of zero. Then, 
we add the two equalities by easily inferring 



2 \\Xn 



H + ^l|VXn(<)||lr + J\\V X n{s)\\ 2 H ds + J (£ n ( S ), X n(s)} ds 



+ / (u n - log u n )(t)dx + / / (1 + Xn)|Vlogu n | dxds + j n / / (1 + X„)(w„- 1) dxds 
Jn JoJn JoJt 

< 77llXOn||ff + ^-W^XOuWh + / (>0r> - log U 0n ) dx 

2 2n J n 

+ [[ (h(9 n ) - \ogu n )x n dxds - [ [ u n Vx n ■ dxds. 



JoJn JoJn u n 

Owing to (13. IS j) and the monotonicity of dyyJ, we have (£ n ,Xn) > a.e. i n (0, T 1 ): indeed, note 
that the value corresponds to the minimum of J (cf. ( 12. 8p ) whence G dyyJ(0). Moreover, since 

< Xn < 1 as well as u n > a.e. in Q (and in E), (4-1) 

and 

r — logr > — (r + | logr|) , r > , (4-2) 
3 

it follows from ([L9]) and (|HTJ) that 

/ / (h(6n) -log u n ) X ndxds < c h + / / | \ogu n \ dxds < c h + 3 / / (u n - log u n ) dxds, 
JoJn JoJn JoJn 

while we remark that 



t r w„. rt 

n 



Vu rt 

u n Vxn ■ — y ^ds = — I I Vxn • V log u n dxds 
u n JoJn 

-oil \^Xn\ 2 dxds + )- 1 1 |Vlogu n | 2 dxds. 
2 JoJn 2 JqJq 



We also recall (I3.5p - (l3.6p and finally apply the Gronwall lemma to deduce that 

\\Xn\\L°°(0,T;H)nL 2 (0,T;V) H 7= || Xn || L°° (0,T; V) + || u n || L°° (O^L 1 (fi)) 



+ || log^n|U°°(0,T;L 1 (n))ni 2 (0,T;V) + 7n||w n ||z,i(£) < C- 



(4.3) 



Second a priori estimate. Test ( 13.1 7ft by <9 t Xn: as the right hand side h(9 n ) — logu n of ( 13 . 1 71) is 
bounded in L 2 (Q) thanks to ( II. 9p and (14. 3p . by a standard computation we deduce that 

||Xn||ir L (0,T;if)nL°°(0,T;V r ) H 7= II ^tXn || L 2 (0,T;V) + ll^(Xn) lU^T^n)) < C- (4-4) 



Third a priori estimate. Rewrite ( 13.171) as 



-A(d t Xn) + A Xn + £n = fn (4-5) 

n 
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with f n := h(6 n ) — logw n — /i d t Xn already bounded in L 2 (0, T; H). Then we can formally test (14. 5 p 
by Axn (such estimate is performed rigorously in [HI formulas (173)-(175)]) and recover, with the 
help of well-known elliptic regularity results, 



IIXn||L 2 (0,T;iy) + II ^Xn || L°° (0,T;H) < C 

provided (13.61) is used to control n~ 1/ ' 2 ||^4xon||_H"- 

Now, thanks to (14.4ft and (14 . 6 j) . it turns out that Cn — fn — MptXn)l n — ^Xn fulfills 

||£n||z, 2 (0,T;V") < C. 

Fourth a priori estimate. We test ( I3.19P by p n and arrive at 



(4.6) 



(4.7) 



\Pn\ 



In 



dxds+~\\V(l* Pn )(t)\\ 2 H +^ / \(1 * Pn )(t)\* dx 



< ||M „||y'||(l *Pn)(t)\\ V < \\u n\\v> / \\p n (s) \\ H ds + \\u Qn \\ V > || V(l * p n ) (t) \\ H 

Jo 

whence, as (1 + Xn) 1 > 1/2 a.e. in Q, using (13. 5j) . Holder's and Young's inequalities we obtain 

'l*Pn)(t)\ 2 dx<C, 



1 

4 



lbn||^ + i||V(l*p n ,)(t)||^+ n 







1 

4" ' v " ™-"- ' 2 
thus without exploiting the boundary contribution on the left hand side. 
Then, recalling also (I3.20p we deduce that 

||Wn||L 2 (0,T;ff) + ||1 * Pn \\ H 1 (0,T;H)nL°°{0,T;V) + \fln~ P * Pn || L°° (0,T;L 2 (r)) < C. (4.8) 

Fifth a priori estimate. We test (13. 22f) by v(l + Un)" 1 where v is an arbitrary function in 
W 1,N+1 (Q) C L°°(Q) and integrate only over Q. Note that, owing to (13. 13[) and the positivity of u n , 
the function v(l + w n ) _1 actually belongs to V. 

Before proceeding, we also remind ( 14. 3 p and immediately deduce that 

|| log(l + u n )\\ L oo (0)T . L i mnL 2 i0:T . v) < C (4.9) 

from the analogous estimate on logti n . 

Then, with the help of (13T22|) and flU) we obtain 



dt loe(l + u„) v dx 



V(U„(1 + Xn)) ■ V 



1 + U r 



dx + 7„ 



< 



< 



(1 + Xn)vdx 

I ((1 + Xn)Vu n + U n VXn) • ( VU - ^ , 2 VM W | rfx 



r 1 + u r 

V 



+ 27„ / \v\ dx 



|Vu I 

(1 + Xn)f- — |Vu| + (l + Xn)M 



|VM r 



;i+w n ) 2 



r/.r 



+ 



— |VXn| |Vu| + — H |VXn 



1 + u 



1 + it„ 



1 + u r 



+ C||u||x,i(r)- 
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Consequently 

\(d t \og(l + u n ),v)\ 

< 2||Vlog(l + Ti n )|| H ||Vv||ff + 2||Vlog(l + « n )||yT;||£cc ( n) 

+ || Vxn||ir||Vu|| Jf + || Vxnlkl|Vlog(l + u n ) \\h\\v\\ l ^{q) + C||tj|| x i (r) 

< C (1 + ||Vlog(l + u n )f H + WVxnfn) lk||^+i ( n) 
a.e. in (0,T). Therefore, in view of (14. 9[) and (14.41) we infer that 

||^log(l + u n )\\ L i {0tT . {w i.N+i mi) < C. (4.10) 

Sixth a priori estimate. We follow the argument devised by Boccardo and Gallouet in [2] and 
introduce the functions r k : R — > R, k G N, defined by 

-1 if r < -{k + 1) 

r + k if - (k+ 1) < r < -k 

if -k<r <k 
r — k if k < r < k + 1 

1 if k + 1 < r 

Note that r k > 0, |r£(r)| < 1, and < r' k ' < 1 a.e. in R. In view of (137L61) and we have that 

Ve n , = (1 + xJnh"(9 n ))V6 n - (h(9 n ) - 9 n h'(9 n ))V X n 

whence 



r fc (0) = 0, r^(r)=<^ 



V0, 



Ve n + (fc(0 n ) - 9 n h'(e n ))V X n 



(1+Xn9nh"(9 n )) 

the denominator being uniformly positive (in fact, bounded from below by due to (11.101) . 

Then, taking r' k (e n ) as test function in (I3.15p . and integrating in time, we have that 

/ T k (e n (t))dx+ / / Tk(e n )Ve n -V9 n dxds 
Jq JoJn 

< / T k (e 0n )dx+ / / T k (e n ) (fi(d t Xn) 2 ~ h{9 n )d t Xn) dxds. (4.11) 
Jn JoJn 

With the help of Young's inequality and fjl.9)) — fll.lOp . we infer that 



/ / T k \e n )X7e n -X79 n dxds 
JoJn 



> r.(e n 



\Ve n \ 



>oJn 

i r* 



l/ (l + Xn9 n h"(9 n )) 
|Ve„| 2 



dxds 



oJn " y " (1 + Xn9 n h"(9 n )) 
1 f*f {h{9 n )-9 n h'{9 n )y< 



(1 + Xn9 n h"(9 n )) 2jJ n (1 + Xn9 n h"(9 n )) 

> ^JJj' k \e n )\Ve n \ 2 dxds - C j\\Vxn\\ 2 H ds. 
On the other hand, we observe that 

fl r' k (e n ) (/i (dtXn) 2 - h(9 n )d tXn ) dxds < C f(l + \\d t x n \\ 2 H ) ds. 
JoJn Jo 



| Vxn| 2 dxds 
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Hence, since Tjfc(e 0n ) < | e n, j and ( 13. 5 p holds, it follows from ( 14. lip and (14. 4p that 

/ T k (e n (t))dx+ [ [ r' k \e n )\Ve n \ 2 dxds<C. (4.12) 
Jn JoJn 

Since |r| < 1 + 7"o(r) for r G 1, it follows in particular from ( 14. 12ft with k = that 

||eri.|U°°(o,T;ii(n)) < SU P / (1 + T (e n (t))) dx < C. (4.13) 
te(o,T) Jn 

Now, if we define D k := {(x,t) e Q : k < \e n (x,t)\ < k + 1}, (T4TT21) yields 

| Ve„| 2 cfeds < C for all keN. 
Taking 1 < q < 2, we infer from the Holder inequality that 



We n \ q dxds = ^2 \We n \ q dxds 



\{2-q)/2 



< El// |Ve.| 2 ^ S \D k f~^<CY,\D, 

k=o \^ D k / fc=0 

Moreover, let g < (iV + l)/(iV + 2) and r = (iV + l)q/N. As k < \e n \ a.e. in D k , we have 









I / |en| 



and, using (I4.13P to control the integral on D , we obtain 

y J |Ve n | 9 dzcfe < C f 1 + J] ^^72 [J J ^ 



fe=i 

k ff/ 2 / oo „ „ \(2-3)/2> 



X] hr(2-q)/q E 



<C|1+I> — - 1 [> // |e n | r <£rds 



k r ( 2 -i)/i 

K k=l / \fe=l 



< C I 1 + ||ej J 2 ^ 72 ( ^ 



l 9/2 N 



k r ^l)ll 

k=l 



r (o — q) 

Let us point out that the series above converges if > 1, which follows from our choices of r 

q 

and q. Therefore, it results that 

\Ve n \ q dxds <c(l+ ||eJ2l 2 (3f 2 ) . (4.14) 
Now, by Holder's inequality we deduce that 

|| e n(*)||L'-(f2) — ll e n(^)|lLi(Q)||en(^)||^5* ( - f2 ^ 

with q* = Nq/(N -q),a + (l- a)q* = r, i.e., 1 - a = (N - q)/N = q/q* if N = 2,3 and q* = oo, 
a — 1, q*(l — a) — r — 1 if N — 1. We use the bound (14.13D to derive 
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and, by Sobolev's embeddings, (I4.13p . and (I4.14p . we have 

||e„ir L , (Q) < C ^ ||e n (t)||^ (Q) dt < C (l + J \\e n {t) - (e re (t), lW^dtj 

- C { 1 + Io l|Ve " (t)l1 W*) " ° v 1 + ll e "ll?'(Q) /2 ) 

from which a bound on ||e n ||£r(Q) follows since 2 — q < 2. Recalling (14. 14ft and the fact that q varies 
in [1, (N + 2)/(iV + 1)) and consequently r in [1, (N + 2)/N), we have thus established that 

N + 2 N + 2 

for all 1 < r < — , 1 < q < ^ + - there holds ||e n || L r ( Q )ni , (0>T;H/ i, 9(c)) < C, (4.15) 

where the constant C depends on r and q as well. 

In view of ( II .9p . (jl.lOp . (]3.16j) . and (14.3)) . it is clear that n satisfies the same estimates: 

for all 1 < r < N ^ 2 , 1 < g < there holds ||6» n || L r (Q)nL9(0iT;H/ i, 9(n)) < C. (4.16) 



Thus, by comparison in (I3.15P and owing to the embedding of L 1 (f2) in (W l ' q> (fi))' , it turns out that 

q N + 2 

— >iV + 2 for all 1< g < — — . 
q-1 ~ N + 1 



ll 9 t e n|Li( ,T;(Wi.«'(n))') < C With V' = Z ~ > N + 2 fOT a11 1 < <? < ( 4 - 17 ) 



This concludes the derivation of uniform estimates for our problem. 



5. Passage to the limit as n ->• oo 

Taking advantage of the bounds in (Q]) - (jPj) . ( 14T6|) (|4~8|) and (j4.15D - fl4.17p . we are in a position 
to pass to the limit as n — > oo by weak and weak star compactness arguments. Hence, possibly 
taking a subsequence (still labelled by n), it turns out that there are functions e, 9, x, £, u , P such 
that 

e n -> e and 9 n -> 9 weakly in U\Q) n L q (0, T; W 1 ' 9 ^)), 

tor all 1 < r < and 1 < o < , 5.1) 

N H N + 1 ' v ; 

weakly star in H\0,T; H) nL°°(0,T;V) C]L 2 (0,T;W) H L°°(Q), (5.2) 

-9 t Xn -> strongly in L 2 (0, T; V% (5.3) 
n 

£ n -»• £ weakly in L 2 (0, T; V"), (5.4) 

u n — > u weakly in L 2 (Q), (5.5) 

l*p n ^l*p weakly star in H l (®,T-H) n L°°(0,T;V). (5.6) 

We aim to prove that (e, 0, x, 6, is a weak solution to problem (jl.ip - (jl.8p fulfilling the conditions 
stated in Theorem 12.21 

Thanks to fl53J and to the boundedness ( jTTD of (<9 t e n ) in L^T; (W 1,9 '(fi))'), it follows from 
the compactness result in [161 Sect. 8, Cor. 4] that 
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e n — > e strongly in L q {Q) for all 1 < q < — ^ ^ , (5.7) 

and this convergence actually takes place in L q (Q) for all q G [1, (N + 2)/N) by interpolation with 
(15.11) . In addition, by the Ascoli theorem and the Aubin-Lions lemma [Til p. 58] we also have 

Xn -+ X strongly in C°([0, T];H)n L 2 (0, T; V), (5.8) 
l*p n ^l*p strongly in C°([0,T]; H). (5.9) 

Now, (15.7p ~ (I5.8p and the Lipschitz continuity of ip^ 1 with respect to both variables (cf. (12. 5ft ) imply 
the strong convergence of 6 n = ip~ l {e n ,Xn) to ip^i^e.x) i n L q (Q) for q G [1, (N + 2)/N); hence, (15.11) 
and the uniqueness of the (weak) limit entail that 

9 n = ip-\e n , xn) ^ 9 = i'^(e, X ) strongly in L q (Q) for 1 < q < — (5.10) 

and (I2.19P follows. By virtue of ( I5.5p -( 15TB|) . (13. ip . and (13. 5p . passing to the limit in (I3.19P we recover 
(I2.22p for u and 1 *p. We point out that equality (I2.23P is a consequence of (I3.20p since {u n {l + Xn)) 
weakly converges to u(l + x) i n L l {Q) by (15. 5 p and (15. 8p . 

In view of ( 14.9^ ( 147X0]) . we use once more [HI Sect. 8, Cor. 4] to deduce that, possibly taking a 
subsequence, 

log(l + u n ) — )• ( weakly in L 2 (0,T;V), strongly in L 2 (Q) and a.e. in Q (5.11) 

for some ( G L 2 (0, T; V). Hence, it turns out that 

u n = exp(log(l + u n )) — 1 — > — 1 a.e. in Q, 

which, combined with (I5.5p . yields u = — 1 and 

u n — y u a.e. in Q and strongly in L S (Q) for all 1 < s < 2 (5-12) 

(see, e.g., [HI p. 12]). Observe next that \ogu is well defined a.e. and integrable: indeed, by the 
Fatou lemma one sees that 

logw| dxds = liminf | logw n | dxds < liminf / / | logw n | dxds < C. 



Then, f l5TT2|) and (0~J) entail 

logn n log u weakly in L 8/3 {Q) n L 2 (0,T;V) 

(5.13) 

and strongly in L S (Q) for all 1 < s < 8/3, 

where 8/3 is chosen in a way that L°°(0, T; L l {Vt)) n L 2 (0,T;V) C L 8 / 3 (Q) for N G {1,2,3}. We 
also recall (I5.10p and the Lipschitz continuity and boundedness of h, which imply 

h{6 n ) -> h{6) strongly in L S (Q) for all 1 < s < oo. (5.14) 

Then, using (I5.2p -( 15T4"P we can pass to the limit in (I3.17P and obtain (I2.20p . From (I3.18P and the 
inequality 

limsup/" (Ut),Xn(t))dt< I (£(t), X (t))dt, 



which is a consequence of (15. 4p and (15. 8p . the inclusion (I2.2ip results from [TJ Lemma 1.3, p. 42] for 
dyyiJ induces a maximal monotone operator also from L 2 (0, T; V) to L 2 (0, T; V). Note in particular 
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that x satisfies (12.271) (look at (I5.2p and use interpolation for the continuity from [0, T] to V) and £ 
is in L 2 (0,T; H), that is (12.281) . due to a comparison of terms in ( 12.201) (cf. also Remark 12.31) . 

Now, it is our intention to show that 

dtXn^dtX strongly in L 2 (Q), (5.15) 

in order to be able to pass to the limit in (13. 15|) . Test ( 13 . 1 T[) by d t Xn an d observe that 

fj, limsup // \d t Xn\ 2 dxds 

n-»oo JjQ 

= limsup (-- // \V(d an )\ 2 dxds -\\\V Xn (T)\\ 2 H - J(Xn(T)) 

n->oo \ Tl J J q Z 



^llVxonlln + J (Xon) + J J (h{9 n ) - log u n )d t Xn dxds 



<-liminf (l\\V X n(T)\\ 2 H + J(Xn(T)) 

- ^W^XouWh - J(Xon) ~ J J (h(9 n ) - log u n )d t Xn dxds 

< - Ql|Vx(r)|||r + J(X(T)) - ^llVxolll - J(Xo) ~ j j W) - \ogu)d tX dxds^ (5.16) 

by the weak lower semicontinuity of the semi-norm in V and the convex function J, as Xn{T) — > x{T) 
weakly in V due to Q and and fET5j» . (l3T7j) . (15TT3|) and (l5TTij) as well. Now, we can exploit 

(I2.20p and conclude that the last line of (I5.16P is nothing but ||^x|Iz,2(q), as one can see by testing 
(I2.20p by dtX and taking advantage of (I2.2ip and [7J Lemma 3.3, p. 73]. Hence, we have that 

limsup || dtXn || l2 (q) < \\dtx\\ 2 L 2 {Q) 

n— >oo 

which, along with the weak convergence (15.21) of {dtXn) to dtX i n L 2 (Q), ensures (I5.15p . 

At this point, we can take the limit in ( 13~T5|) . Thanks to (j3T5]l . f l5Tj) . (QUI) . (15TT4D - (15TT5D it 
turns out that 

e„(t) = e „ - / (A0 n + h(6 n )d t Xn ~ ^{dtXnf) (s) ds 



converges in a suitable way to 

e(t) = e - / (A0 + /i(#)<9 a - /i (9a) 2 ) (a) <fc, 

whence the initial condition (12.241) for e and (12.181) . Note in particular that e is continuous from 
[0,T] to (W 1 > q ' {Q))' and satisfies (T2TT3|) and (f2~26]) for all r, g fulfilling 1 < r < (TV + 2)/N, 1 < 
g < (N + 2)/(N + 1). In fact, the latter comes as a consequence of the bounds (I4.15P coupled 
with the convergences (15.11) . About the L°°(0,T; L 1 (fi))-regularity of e, (15.71) implies that e n — > e in 
^(0,7; L 1 ^)), whence 

||e n (£)||.Li(fi) — > ||e(t)||z,i(n) for a.e. £g(0,T), at least for a subsequence. 
Then, recalling (I4.13P we infer that 

IKOIU 1 ^) - SU P 1 1 ^ 1 1 (o.rji 1 (ri)) < C, 
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and consequently e G L°°(0, T; The same property can be deduced for u and logw (cf. (I5.12p 

and (15 . 131) ) so that (12.291) holds as well. The proof of Theorem 12.21 is then complete. 



6. Complementary results 

This section is devoted to the proofs of the other results stated in Section [2j We start with 
Theorem 12.41 ensuring the positivity of the variable 9. 

Proof of Theorem 12.41 Let us go back to our approximating problem ( 13.15l) -( !3.21l) . recalling in 
particular Proposition 13.11 By (SJ Theorem 2.6] we have 

log^GL oo (0,T;L 1 (fi))nL 2 (0,T;\/), % G L 2 (Q) (6.1) 
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whenever log^Qn G L 1 (f2), besides the other requirements ( I3.2p -( 13~4"|) on the initial data. Then, here 
we have not only to guarantee that log#on G L 1 (f2) for all n G N, but we also need the uniform 
boundedness property 

||log0on|k(n)<C (6.2) 
along with (I3.5p -( l3~6|) to prove (I2.32p for the limiting functions. 

By (12.41) . (12. lip , and (I2.3ip we know that both 8q and log#o are m L 1 ^). Then we could, for 
instance, introduce r] := y/9o, observe that r] G H and logr/o = 2 log^o £ ^(Q), an d consequently 
approximate r] by taking the solution r] 0n to (cf. ( 13.91) ) 

t] 0n + -Ar] 0n = 7] + - in H, (6.3) 
n n 

then setting 9 0n := ?7 2 n . By this choice, it turns out that T] 0n G W (whence 9 0n G V) and rj 0n > 1/n 
a.e. in Q. Moreover, r] 0n — > r] in H so that 8 Qn — > O in L 1 (fi) as n — > 00. Then the related 
sequence e 0n := ip(9 0n ,Xon) satisfies (I3.5p i . i.e., converges to e in L 1 (fi) thanks to the properties of 
if) (see (I2.4p -( l2~5|) ). Furthermore, if we test (16. 3 p by (1 — l/?7o n ), use (I4.2p . and proceed as in the 
computation following (I3.9p . we arrive at 

\ I |log^on|^<o / (\von\ + \\ogr] 0n \) dx 

< / (VOn ~ log?7on) dx < / (^ ~ log T]o) dx + | fi| , 

Jn Jn 

whence (I6.2p is satisfied. 

Now, we can deal with the approximating problem as in Proposition 13.11 and, in particular, consider 
equation (I3.24[) which is an alternative formulation of (I3.15p . We formally test (I3.24p by —l/8 n and 
integrate by parts: we easily infer that 

- log 6 n )(t) dx + 1 1 — TTj- dxds + fi I f ■ t -^ L dxds 
n JoJn JoJn 

< / (- log6» „) dx + / (h'(6 n )xn)(t) dx - / h'(6 0n )xondx. 
Jn Jn Jn 
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Then we add the term / 9 n (t) dx to both sides of the above inequality and recall (11.9j) . (14. II) . ( I4.2j) . 

in 

and (I4.16P to deduce that 

\ [ (K(t)\ + \ log 9 n (t)\) dx + [ [ |V log 6 n \ 2 dxds + fi [ [ ^X^ 2 dxds 
3 Jn JoJn JoJn "n 

< f(\\og9on\ + \h'(9 0n )xon\)dx+ f(9 n + h'(9 n ) X n)(t)dx<C, 
Jn Jn 

whence we find out the estimate 

dtXn 



l°g 9 n 1 1 Loo (o,T;Ll (H))nL 2 (0,T; V) + 



'9, 



< C. (6.4) 

i 2 (Q) 



Next, we can pass to the limit as n — > oo. In view of (15. 10p and (15. 15j) . we have that 9 n — > 9 and 
dtXn — > dtX a -e. in Q, at least for a subsequence. Thus, similarly as in (15. 1 3[) we infer that 9 > a.e. 
in Q and 

\og9 n \og9 weakly in L 8/3 (Q) nL 2 (0,T;V) 

and strongly in L S (Q) for all 1 < s < 8/3, 

— = weakly in L 2 {Q) and strongly in L S (Q) for all 1 < s < 2. 



vC V9 

At this point, it is a standard matter to recover (j2.32j) and Theorem 12.41 is completely proven. □ 

The next step is showing the stability and decay estimate in the case 7 > 0. 

Proof of Theorem [2J3 We first note that for a fixed x fulfilling (I27TI]) . ([2725]) and ( 127271) (as it 
is for our solution), there exists a unique pair (u,p) solving (12 . 2 2 j) (12 . 23 p and such that 

ueC°([0,T];V')nL 2 (0,T;H), (6.5) 
peL 2 (0,T;H) with l*peC°([0,T];V) (6.6) 

for all T > 0. Now, we proceed formally by letting u and p be smooth enough in order to perform 
our computations: for instance, a rigorous proof can be obtained by suitably regularizing uq and the 
fixed x an d then passing to the limit. Hence, we consider the solution ( G C°([0, +00); V) to 

AyC(t) = u(t) in V, for all t > (6.7) 

and perform a formal estimate by testing ( 16. 7p by dt((t). Thus, with the help of the time differentiated 
version of (12722]) . 

d t u{t) + A^t) = 0, te(0,T), 

we easily obtain 

^<AC(t).C(*)> = (ACtW.CW) = <M*U(*)> = -(A, P (t)x(t)) 

= -(A 7 C(t),p(t)) = -(u(t),p(t)). (6.8) 

Owing to ( I2.23H and setting 



^(<):=^l|VC(*)||i + ^ / ICWI 2 ^, $ 2 (t):= ff (l + x)u 2 dxds, 



r Jo-'n 
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from (16. 8p it turns out that the continuous and nonnegative function $ x + $ 2 is constant, or equiva- 
lent ly, 

^P~(t)+ [ ((1 + X )u 2 ) (t) dx = for a.e. te(0,+oo). (6.9) 
Now, since (16 .7p implies 

1 /2 

I|VC(0II^ + 7^ICWI 2 ^< \Ht)\\v>\\((t)\\v < c \\u(t)\\ v , (||vc(t)||^ + 7^ICWI 2 ^ 

and, recalling that 

\\u(t)\\ 2 v , < C\\u(t)\\ 2 H < C [ ((1 + xK) (t)dx 



for a.e. t G (0, +oo), it is not difficult to see that (I6.9P entails 

^■$i(t)+a*i(t)<0 (6.10) 

for some positive constant a. Then the Gronwall lemma yields < $i(0)e~ at , which means 

exponential decay for both ||C(0llv an d |l M (^)||v since (cf. (16. 7p ) 

= H^rCWllv < C7||C(*)||v- 

Moreover, a comparison in (I2.22p shows that (l*p)(t) converges to A~ 1 («o) in V (still at exponential 
rate) as t — > +oo. Now, by integrating (16. 9 j) with respect to time and using the nonnegativity of x, 
we have 

\\u\\L 2 (o,t;L 2 (Q)) = u 2 dxds < $ 2 (t) < $i(0) < C for all t G (0, +oo), 

which completes the proof of Theorem 12.51 □ 

Finally, we prove the result on steady state solutions for the case 7 = 0. 

Proof of Theorem 12.61 Here, we are looking for steady state solutions of problem f l 1 . 1 1) — ( TXTHT) 
with 7 = 0: let us remind Definition 12.11 and look for sextuples of functions (e, 8, x, u , p) which 
do not depend on time. 

In view of (12.181) . we observe that A9 = implies that 9 is a constant function; the same can be 
concluded for p from (I2.22p . We now look for stationary \ solving (I2.20p . which can be conveniently 
rewritten as 

Ax + co = h{6) - logp (6.11) 
where (cf. (II. 4p ) oj := £ — log(l + x) satisfies 

u) G dl[ 0i i](x) a.e. in f2. (6.12) 

Since the right hand side of (16. lip is constant, two solutions xi? X2 of (16. lip may differ just by a 
constant: in fact, if we test the difference of equations by Xi ~ X2, due to the monotonicity of dlmii 
we have that |V(xi — X2) I — a.e. in Q. Therefore, we can observe that (recall that p = u(l + x)) 

(i) if h(8) — logp > then x — 1 an d oj = h{9) — logp solve (I6.lip -f l6.12p . According to the 
above discussion, any other solution x should be a constant, but no constant different from 1 
is allowed by (I6.12p . Therefore, x — 1 an d u — h(8)—logp is the only solution to (I6.1ip - (l6.12p 
in that case. 
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(ii) if h(9) — logp < then \ = and u = h{9) — logp solve (I6.11l) - (l6.12p . Any other solution 
X should be a constant, but no constant different from is allowed by (16.121) . so that \ = 
and cu = h{6) — logp is the only solution to ( 16.11l) -( |6.12l) in that case. 

(iii) if h{6) — logp = then any constant x £ [0, 1] solves (I6.1ip - (l6.12p with uj = 0. 

Then x turns out to be a constant in all cases and consequently also e = ip(6, x) given by ( 12.191) is a 
constant. This concludes the proof of Theorem 12.61 □ 
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